A One-sided, Highly Relativistic Jet from Cygnus X-3 
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O ■ ABSTRACT 
H ■ 

' Very Long Baseline Array images of the X-ray binary, Cygnus X-3 were obtained 2, 4 and 

Ci • 7 days after the peak of a 10 Jy flare on 4 February 1997. The first two images show a curved 

one-sided ]ei, the third a scatter-broadened disc, presumably at the position of the core. The jet 
curvature changes from the first to the second epoch, which strongly suggests a precessing jet. 
f\ ' The ratio of the fiux density in the approaching to that in the (undetected) receding jet is > 330; 

^ . if this asymmetry is due to Doppler boosting, the implied jet speed is > 0.81c. Precessing jet 

model fits, together with the assumptions that the jet is intrinsically symmetric and was ejected 
during or after the major flare, yield the following constraints: the jet inclination to the line of 
sight must be < 14°; the cone opening angle must be < 12°; and the precession period must be 
> 60 days. 

Subject headings: binaries: close — stars: individual: Cygnus X-3 — radio continuum: stars 

1. Introduction states, Cygnus X-3 is the most luminous X-ray bi- 

nary at radio wavelengths (Waltman et al. 1995). 
Cygnus X-3 is one of the few X-ray binaries j^s companion is thought to be a Wolf-Rayet star 

that is consistently strong in both radio and X- (^^n Kerkwijk et al. 1992), and a 4.8 h period has 

ray emission. In fact, in both quiescent and flare ^^^^ seen at X-ray (e.g.. Mason & Sanford (1979)) 

; and infrared (Mason, Cordova, & White 1986) 

1 JIVE, Postbus 2, 7990 AA, Dwingeloo, the Netherlands . ^, ... . , ' 

T .T,,.^ T.,^ ^ r. energies, ihis period is generally interpreted as 

^present address: NRAO, P.O. Box 0, Socorro, NM , i • , • , n i i • m , 

87801 USA orbital period oi the binary system. Schmutz, 

^deceased, July 29, 2000 Geballe, & Schild (1996) used the velocity shifts 

of lines in the infrared to estimate, assuming the 
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orbital system is seen edge-on, that the compact 
object is a black hole with mass > 7 Mq. 

Giant radio outbursts have been known in 
Cygnus X-3 since 1972 (Gregory et al. 1972). 
Large radio flares are always preceded by a 
quenching of the radio emission (Waltman et al. 
1994). The hard X-ray is usually anti-correlated 
with the radio except during quench-flare peri- 
ods, although the quench period in hard X-ray 
is generally longer than in the radio (McCoUough 
et al. 1999). Similarly, the soft thermal X-ray and 
radio emission is usually correlated except dur- 
ing quench-flare periods (Watanabe et al. 1994). 
In other words, before a large radio/hard X-ray 
flare the soft X-ray is more luminous than usual, 
and during a flare the soft X-ray intensity briefly 
drops. Observations with the VLA (Geldzahler 
et al. 1983), MERLIN (Spencer et al. 1986), and 
VLBI (Schahnski ct al. 1995; Molnar, Rcid, & 
Grindley 1988; Schalinski et al. 1999), during 
or shortly after large radio flares, suggest ejection 
of radio-emitting plasma in the north-south di- 
rection. Spencer et al. (1986), Schalinski et al. 
(1995), and Molnar, Reid, & Grindley (1988), 
assuming a two-sidcid jet, estimated an expansion 
rate of ~ 5 mas/day. This corresponds to an ap- 
parent transverse speed of 0.3c for a distance of 
10 kpc, which we adopt for the remainder of the 
paper (based on H I absorption; see Dickey 1983). 
However note a recent H I absorption observation 
by S. J. Befl BurneU & W. M. Goss (priv. comm) 
which implies a distance of 11.5 ± 1 kpc based 
on i?o = 8.5 kpc; and, a Chandra observation of 
Cygnus X-3's X-ray scattering halo from which 
a geometric distance of 9 ±2 kpc was calculated 
(Predehl et al. 2000). 

Here we present high resolution images of 
Cygnus X-3 taken with the Very Long Baseline 
Array (VLBA)'', following a large radio flare in 
February 1997, showing a one-sided radio jet di- 
rected south from a highly- variable core. §2 de- 
tails the observations and data reduction, while 
§3 describes the images and model fits, and §4 the 
conclusions. 



^The VLBA is a facility of the National Radio Astronomy 
Observatory which is operated by Associated Universities, 
Inc. under cooperative agreement with the National Sci- 
ence Foundation 
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Fig. 1. — Radio light-curves for Cygnus X-3 at the 
time of our VLBA observations. The horizontal 
lines show the times spanned by our VLBA obser- 
vations. Note that the source was optically thin 
during the first epoch, and optically thick towards 
the end of the third; there arc no simultaneous flux 
measurements for the second epoch. The x-axis is 
Modified Julian Date: MJD= JD-2400000.5. 

2. Observations and Data Reduction 

Cygnus X-3 was observed with all 10 anten- 
nas of the VLBA in three 13 hour sessions, two, 
four, and seven days after the source flared to 
above 10 Jy at 15 GHz on 4 February 1997. Fig- 
ure 1 shows the epochs of the VLBA observa- 
tions in the context of the radio light-curves, taken 
from the Green Bank Interferometer (GBI) and 
the Ryle Telescope (to be discussed in a later pa- 
per) . During most of the first epoch (6 February) 
observations the source as a whole was optically 
thin (a ~ —0.6, Su oc v°') and rapidly decaying, 
from 2.1 to 1.75 Jy at 8.4 GHz. Unfortunately 
there were no total flux density measurements dur- 
ing the second epoch (8 February); by the third 
(11 February) Cygnus X-3 had begun a series of 
smaller flares with quite rapid declines, the flux 
density at 15 GHz varying by a factor three or 
more during the VLBA observations. 

A third and potentially more serious complica- 
tion is the motion of the jet during the observa- 
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Fig. 2.- VLBA images of Cygniis X-3, from 6, 8, and 11 Fcbniary 1997. The restoring beam (a Gaussian 
fit to the dirty beam) is shown in the lower left corner of each image, a. 6 February 1997. Contours are 
±20 X 2"/^ mJy/beam, n = 0,1,2,...; the peak in the image is 161mJy/beam, with the northern core 
roughly half that bright. The gray scale spans the range of the image, from —29.4 to 160.8 mJy/beam. The 
rms noise is 7.5mJy/beam, and the restoring beam is 3.98 x 3.00 mas. b. 8 February 1997. Contours are 
±1.25 X 2"/^ mJy/beam, n = 0,1,2,...; the peak in the image is 33.9 mJy/beam. The gray scale spans 
the range of the image, from —2.0 to 33.9 mJy/beam. The rms noise is 0.48 mJy/beam, and the restoring 
beam is 4.66 x 4.13 mas. As discussed in the text, the time- variable core has been removed and replaced 
by a Gaussian with its mean flux density, c. 11 February 1997. Contours are ±60 X 2"/2 mJy/beam, 
n = 0, 1, 2, . . .; the peak in the image is 322.4 mJy/beam. The gray scale spans the range of the image, from 
—73.5 to 322.4 mJy/beam. The rms noise is 20 mJy/beam, and the restoring beam is 4.29 x 3.81 mas. Unlike 
the 8 February image, no attempt has been made to remove the time-variability of the core in this image. 



tions. As discussed below, the proper motion of 
the jet is > 20 mas/day, corresponding to mov- 
ing by two beams or more during our ~ 1/2 day 
observations. 

The VLBA observations were carried out with 
a 16 MHz bandwidth centered on 15.365 GHz, us- 
ing dual polarization and two-bit sampling. The 
data were correlated using the VLBA correla- 
tor in Socorro, NM, USA, and amplitude cali- 
bration and fringe fitting were performed using 
NRAO's Astronomical Image Processing System 
(AIPS). Scans on Cygnus X-3 were "sandwiched" 
between scans on the calibrator J2025±334, 7° 



away, with a three-minute cycle time, allowing the 
use of antenna amplitude and phase gain solutions 
from that source as a first estimate for those on 
Cygnus X-3. 

There were several additional complications to 
the data reduction. First, J2025+334 was too 
far away from Cygnus X-3 to track the latter's 
phase variations accurately; although the position 
of the northern-most component (presumably the 
core) was stable to within a few milliarcseconds, 
self-calibration was required to bring out the full 
structure of the jet. Second, both Cygnus X-3 and 
J2025+334 are scatter-broadened, as indicated by 



the fall-off in the interferometric amplitudes on the 
longest baselines even when the core completely 
dominates the image (as in the third epoch). This 
may of course be intrinsic source structure, but 
is consistent with the scattering sizes previously 
measured (see §3.1). In general we concentrate 
on the more extended jet emission, and the im- 
ages of all three epochs, shown in Figure 2, are 
made with a 30MA full-width at half-maximum 
(FWHM) Gaussian taper. The corresponding syn- 
thesized beams are roughly 3-5 mas FWHM. 

Simulations showed that proper motions of this 
order would have only minor effects on the im- 
ages, and would not affect any of our major con- 
clusions. In particular, we modeled the source as 
a stationary core with a straight jet 60 mas in 
length beginning 10 mas away, oriented along a 
position angle of 175°. The total flux density in 
the jet was taken to be about 10 times that of 
the core, and fell off as where r is the dis- 

tance along the jet; this roughly matches the im- 
age from the first epoch (Figure 2a). The simu- 
lated jet was taken to move outward as a whole 
along the same position angle at 20 mas/day. The 
core, and the jet perpendicular to the ejection axis, 
were modeled as unresolved, but convolved with 
a Gaussian of 2 mas FWHM to simulate the ef- 
fect of interstellar scattering. The Fourier trans- 
form of this model was then projected onto the 
observed baselines in the uv-plane, and the entire 
13 hour data set imaged and deconvolved follow- 
ing the same procedure as used for the observed 
data. The resulting image showed a slightly elon- 
gated version of the model as it appeared at the 
mid-point of the observations, with the elongation 
corresponding to the proper motion of the jet over 
13 hours; this seems intuitively reasonable, and 
shows that proper motion does not significantly 
affect our conclusions. Following individual radio 
components between epochs is not advisable, but 
the apparent curvature and rough extent are real. 

The February 6 data were further compromised 
by poor a priori amplitude calibration of several 
antennae due to snow over most of the south- 
western United States. This epoch required corre- 
spondingly more drastic self-calibration and care- 
ful imaging, and the resulting map (Figure 2a) is 
still not nearly as good as those from the second 
epoch. 

The 8 and 11 February data were taken in 



much better weather, and for these epochs the ma- 
jor complication is the intrinsic variability of the 
source. This produces artifacts which dominate 
the noise in the maps, and makes imaging even the 
core itself quite difficult. During the second epoch 
(8 February) most of the variability, as well as the 
total flux, was concentrated in the core. We there- 
fore created a core-only data set by subtracting a 
first (crude) model of the jet, self-calibrated that 
data set in hour-long segments, transferred the an- 
tenna gains thus derived to the original data, and 
subtracted CLEAN-component models (Hogbom 
1974) of the time-variable core. This produced a 
jet-only data set phase-referenced to the core. Fig- 
ure 2b shows the image produced from this data 
set with one round of phase self-calibration applied 
and a Gaussian model of the core with an average 
flux added. The improved calibration reduced the 
rms noise in the image by a factor of 5, showing 
that the jet extends an additional 90 mas beyond 
what could be seen in the maps which were phase- 
referenced solely to J2025-I-334. 

Unfortunately this method did not improve the 
6 and 11 February data sets. For 6 February this 
could be because of the poor initial amplitude cal- 
ibration, or because the flux variability was pri- 
marily in the jet rather than in the core; most 
likely, both effects contribute. The 11 February 
map (Figure 2c) by contrast is entirely dominated 
by the core, but the above procedure while reduc- 
ing the noise level did not reveal any convincing 
jet, to a limit of 10-20 mJy/beam in a 5 mas beam. 
This noise level is far higher than that in the 8 
February map (Figure 2b), presumably because 
the core is much stronger and its flux variations 
far greater (a factor 3-4 during these observations; 
see Figure 1). With this noise level one could not 
expect to see any emission corresponding to the 8 
February jet, even if that emission had not faded 
at all. 

3. Discussion 

3.1. Is the Northern Component the Core? 

One of the primary conclusions of this work 
is that the jet in Cygnus X-3 is one-sided, mak- 
ing this the first severely asymmetric Galactic 
relativistic source. This conclusion rests on the 
premise that the northernmost component is as- 
sociated with the core. The evidence for this 
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is quite compelling. Although phase-referencing 
was not entirely successful, the maps before self- 
calibration did show that the brightest (northern- 
most) component is stationary to within 3 mas. 
Further, the strong variations seen with the Ryle 
telescope on 11 February are mirrored exactly in 
the amplitudes measured on the shortest VLBA 
baselines, consistent with the corresponding im- 
age (Figure 3), which shows no evidence for any 
but the compact northern component. This vari- 
ability, on timescales of minutes to hours, is most 
easily explained by assuming this component is in- 
deed the central core, with the variations caused 
by short-timescale ejection and decay of unre- 
solved jet segments. Finally, simultaneous mea- 
surements with the GBI towards the end of the 
VLBA observations show that the source was op- 
tically thick (a ~ 0), a characteristic of emission 
on very small scales; the fact that the emission at 
this time came solely from the northern compo- 
nent is another argument that that component is 
indeed the core. 

One possible counter-argument is that the 
northern component is clearly resolved, as seen 
both in images and more directly in the uv- 
plane, where the visibility amplitude drops off 
rapidly with radius. This might be either in- 
trinsic structure or broadening due to interstel- 
lar scattering; Cygnus X-3 is one of the most 
scatter-broadened sources known (e.g., Wilkinson, 
Narayan, & Spencer 1994; hereafter WNS), and so 
might be affected even at these high frequencies. 
To check this quantitatively, we used K. Desai's 
private AIPS task OMFIT to fit half-hour seg- 
ments of the 11 February uv-data directly, both 
with and without simultaneous short time-scale 
(10-120s) phase self-calibration. Segmenting the 
data was required to disentangle the source size 
from the extreme flux variability; fitting in the 
uv-plane is preferable to more normal image-plane 
fitting, both because those data are the observed 
quantities, and because the poor uv-coverage in 
such short periods makes mapping and deconvo- 
lution difiicult. In any case the results from OM- 
FIT were consistent with more standard fitting in 
the map plane, and the residuals were reasonable 
both in the map and in the uv plane. The best-fit 
Gaussian has a FWHM (geometric mean of major 
and minor axes) of 1.7 ± 0.3 mas, an axis ratio of 
0.8 ± 0.25, and is elongated along a position angle 



Cyg X-3: Ryle vs. VLBA (BH039C, 10-15MA) 
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Fig. 3. — Comparison of Ryle Telescope and 
VLBA flux densities at 15 GHz, for 11 Febru- 
ary 1997. The solid circles show the Ryle data, 
the open squares (RCP) and triangles (LCP) the 
average of the VLBA baselines between 10 and 
15 MA. The solid line shows the Ryle flux den- 
sities divided by 1.6; the agreement is excellent, 
apart from this constant factor. This multiplica- 
tive offset is due partly to scatter-broadcning (for 
any reasonable image the flux density increases as 
the baseline gets shorter, and the Ryle data are ef- 
fectively taken at OMA), and partly to calibration 
errors (see text). 

of 50 ± 20° (the error bars reflect both statistical 
errors in the individual fits, and the full scatter 
between the fits for the half- hour segments). This 
and previously published size determinations are 
shown in Figure 4. Our measurement, as well as 
the 22 GHz observation of Molnar, Reid, & Grind- 
ley (1988), is consistent with a v~'^ extrapolation 
from earlier low-frequency (< 5 GHz) measure- 
ments of the scatter broadening, but disagrees 
both with 8.4 GHz observations by Geldzahler, 
Kellermann, & Shaffer (1979; hereafter, GKS) 
and with 15 GHz observations by Newell, Garrett, 
& Spencer (1998; hereafter, NGS). The disagree- 
ment with GKS (who obtain a size of 1.3±0.2 mas) 
might be dismissed on the grounds that their re- 
sult was based on a single short observation with 
poor sensitivity and only four useful baselines; 
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their measurement also falls well below any rea- 
sonable extrapolation of the lower frequency data. 
NGS on the other hand observed for 9 hours with 
the full VLBA, producing a data set quite compa- 
rable to ours, and obtained mean FWHMs of be- 
tween 1 and 1.8 mas at roughly the same frequency 
(15.3 GHz); they also found that the FWHM was 
strongly correlated with the flux density on short 
VLBA baselines. (At the current stage of analysis 
our data do not demand that the source change 
size, and in any case show with high confidence 
that the size at all times lies within the bounds 
noted above.) 

It is difficult to reconcile NGS' lower bound 
on the size (1 mas) with the previous measure- 
ments and any reasonable scattering model, since 
it lies significantly below a v~'^ extrapolation of 
those data. On the other hand, the size (1.8 mas) 
and orientation (position angle ^ 60 70°, based 
on their images) NGS found when the source was 
brighter matches our own. Further, both the el- 
lipticity and position angle derived from our data 
are within the errors identical with those derived 
at 1692 MHz by WNS (1.32 and 62°, respec- 
tively) and at 1665 MHz by Molnar et al. (1995, 
1.31 ± 0.02 and 52?0 ± 1?5). Admittedly most of 
the observations, as well as our own, were taken 
during or shortly after flares, and NGS found their 
larger sizes during short-lived 'mini-flares' to 200- 
300 mJy. But while NGS find a much smaller size 
(^ 1 mas) during quiescence [Su ~ 40-60 mJy at 
15 GHz), Fender et al.'s (1995) three-day MER- 
LIN observations at 1658 MHz during another qui- 
escent period (40-50 mJy at both 18 and 2cm) gave 
a size which was only slightly smaller than ear- 
lier flaring measurements (142 ± 5 mas, vs. 150- 
160 mas from Spencer et al. 1986, WNS, and Mol- 
nar et al. 1995 at 1660-1692 MHz) but with a 
virtually identical axis ratio (1.31 ± 0.08) and po- 
sition angle (62° ± 3°). The Fender et al. (1995) 
result seems fairly solid, given the large amount 
of data (three full tracks, yielding an rms noise of 
38/LtJy) and high resolution (0.15 arcsec). If the 
NGS result is correct their observations must have 
coincided with a very unusual period of lower and 
much more isotropic scattering. 

So we are left with two possibilities: either 
the NGS result for some reason is spurious, and 
the northern component is an unresolved source 
broadened by interstellar scattering; or the NGS 



result is correct, and the northern component we 
see has roughly the same size and orientation as 
they observe during mini-flares. Either way it is 
difficult to imagine that the northern component 
is not closely associated with the central object, 
and it seems most reasonable to assume that it is 
the origin of the extensive jet ejected during the 
February radio flare. All the evidence therefore 
indicates that our images of Cygnus X-3 indeed 
show a one-sided jet source. 

Cygnus X-3 Scattering 




Frequency [GHz] 

Fig. 4. — The apparent size of Cygnus X-3 

hmction of frequency. Sizes are geomet- 
ric means of the minor and major axis FWHM 
of elliptical Gaussian fits, or of circular Gaus- 
sians if only those were fit. The solid line rep- 
resents a weighted least-squares fit to all the 
data, e = 448 {v/l GHz)"^ °^ mas; the dashed 
line gives the fit using only data below 6 GHz, 
= 453 (i^/l GHz)^^ °^ mas. Data are taken from 
Anderson et al. (1972); Fender et al. (1995); 
Gcldzahlcr et al. (1983); GKS; Molnar, Rcid, & 
Grindley (1988); Molnar et al. (1995); NGS; 
Schalinski et al. (1995); Spencer et al. (1986); 
WNS; and this work. Two values (represented by 
crosses) are given for NGS at 15 GHz, representing 
the quiescent (minimum) and minor flare (maxi- 
mum) states. The maximum NGS value is overlaid 
by our measurement, which is indicated by a solid 
circle. 
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3.2. The Jet 

Accepting the northern component as the core, 
the jet, shown in Figure 2, is at least 50 mas long 
on 6 February, and 120 mas long by 8 February. 
These are lower limits on the true length of the 
jet, as it fades into the noise towards the south. 
In both maps the jet is curved, and that curvature 
changes between the two epochs. 

What can we learn from these images? Most 
obviously, the ejection must be fairly continu- 
ous, a stream rather than one or a few 'blobs.' 
This is reminiscent of the jets seen in SS433 (Ver- 
meulen et al. 1987) and GRO J1655-40 (Hjellming 
& Rupen 1995), but unlike the easily-separable 
components seen in GRS 1915+105 (Mirabel & 
Rodriguez 1994; Fender et al. 1999a). Unfor- 
tunately this, combined with the poor quality of 
the 6 February image, makes it impossible to mea- 
sure the proper motions of individual components 
directly. 

Based on the radio light curves (Figure 1) 
it seems likely that the jet was ejected at the 
time of the large radio flare, beginning roughly at 
MJD 50482. liO.l (MJD - Modified Juhan Date: 
JD-2400000.5). This is not absolutely conclusive 
— the amount of flux in the most extended struc- 
tures is not so great that it would necessarily have 
dominated the light curve even a week or more be- 
fore the flare, if it did not decay over time. But 
the jet emission has clearly declined significantly 
from the first to the second VLBA images (Fig- 
ures 2a and 2b), and this together with the quite 
sharp decline seen in most of the smaller radio 
flares strongly suggests that the emission we see 
originated in the main radio flare on MJD 50482.1. 

Table 1 lists the epochs of observation, inferred 
ages, proper motions and apparent jet speeds. 
These proper motions (14 — 21 mas/day) are sig- 
nificantly higher than those previously cited for 
this source, which range from 8.4 (Schalinski 
et al. 1995, corrected for a one-sided jet) to 
10 ±2 (Geldzahlcr ct al. 1983) to 4.6-18 mas/day 
(Spencer et al. 1986, with the range correspond- 
ing to the uncertainty in the ejection date). Given 
the limitations of those earlier observations, which 
could easily have missed the sort of low-level ex- 
tended structure on which we base our higher 
estimate, there is no evidence that the proper mo- 
tion has changed with time. To the contrary, the 



agreement between the position angles measured 
in those earlier observations (in all cases almost 
directly north-south) and our own suggests that 
the ejection axis at least has remained remarkably 
stable since at least 1982. 

One possible inconsistency is that the uv-plane 
fits to the data from the third epoch do not suggest 
such large source motions. Probably the expla- 
nation is that the emission associated with short 
timescale flares is intrinsically different from the 
longer-lived emission characteristic of larger flares 
like that of 4 February; for instance, it might be 
that the emission associated with the short flares 
decays too rapidly to be seen outside the scattering 
disk. For the remainder of this section we assume 
that the third epoch provides no relevant informa- 
tion on the motions of the jet seen in the flrst two 
epochs. 

3.2.1. Physical Parameters 

Our images give three basic observables useful 
for deriving the intrinsic properties of the jet: the 
apparent proper motion; a limit on the ratio of 
the brightness of the approaching to that of the 
receding jet, if the system is intrinsically symmet- 
ric; and the appearance of the jet itself, in partic- 
ular its curvature and evolution with time. The 
flrst two are susceptible to direct mathematical 
analysis; the last requires more heuristic modeling, 
which is deferred to the next subsection. Note that 
we implicitly assume that the observed proper mo- 
tion is related to a physical velocity, rather than a 
group speed or some optical illusion involving e.g. 
scattering screens or the simultaneous 'lighting up' 
of well-separated hot spots. 

Under this assumption, according to special rel- 
ativity (e.g., Hughes & Miller 1991), a jet moving 
at an intrinsic speed fic at an angle i to the line-of- 
sight will be observed to have an apparent trans- 
verse motion /3appC, where 

/? sin i 
1 — p cos I 

The minimum /3 for the approching component, 
for a given /3app is 

r, _ Pa,pp 
Pmin — J 
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Table 1 
Observations 



Epoch 

MJD 


date 


time since start of flare 
(days) 


proper motion 
(mas/day) 






50485.46-50486.00 
50487.46-50488.05 
50490.46-50491.04 


(6 Feb.) 
(8 Feb.) 
(11 Feb.) 


(3.4 -3.9) ±0.1 
(5.4- 5.9) ±0.1 
(8.4 - 8.9) ± 0.1 


14.7 ±0.4- 12.8 ±0.3 
21.5 ±0.4- 19.7 ±0.3 


^ (0.85 
> (1.24 


- 0.74) ±0.02 

- 1.14) ±0.02 



Note. — The proper motions are based on lengths of 50 mas and 120 mas for 6 and 8 February, respectively. The apparent 
velocity is based on a distance of 10 kpc. 



corresponding to an inclination of sin y i+p^ 

As discussed above, for our observations /3app ^ 
(1.24-1.14)±0.02, implying > (0.78-0.75). 
Since /3 < 1, the inclination lies between 0° and 



sm 



app 



1 + 0"^ 
' r-'app 



For Cygmis X-3 this implies an inclination be- 
tween 0° and (78 - 83)° ± 1°. The correspond- 
ing Doppler boosting (see below) is unconstrained; 
for inclinations above 59°, the source could ac- 
tually be de-boosted, i.e. appear fainter than it 
would without relativistic effects, as has been seen 
in GRS 1915+105 (Mirabel & Rodriguez 1994; 
Fender et al. 1999a). 

If the source were intrinsically symmetric, the 
ratio of the flux in the approaching to that in 
the receding jet would provide an additional con- 
straint. Since we do not detect a jet on the oppo- 
site side, we have only an upper limit on this ratio: 
-Robs ^ 330, determined by integrating the flux in 
the southern jet and comparing that to the noise 
to the north, integrated over a somewhat smaller 
region to account for the slower (apparent) mo- 
tion of the receding jet. The Doppler factor for a 
relativistic jet is 



D = 



1 — (3 cos i 



The corresponding Doppler boosting of the ob- 
served flux density is (Pearson & Zensus 1987) 

Si, 



where v is the frequency, Sv the flux density. A; is a 
constant (fc = 2 for a continuous jet, k = ?> for an 
unresolved blob) , and a is the spectral index {Su oc 
I/"). The primed quantities refer to the rest frame 
of the object. The predicted flux ratio (between 
the approaching and receding jets) is then 



R 



1 + /3 cos i 
1-/3 cos i 



k—a 



where i is now taken as the inclination of the ap- 
proaching jet. Solving for the intrinsic as a func- 
tion of the observed quantities, we have 



(3cosi 



R-- 



1 



R~- 



1 



Since both /? and cosi can be at most 1.0, R> 1, 
and (k — a) > for synchrotron emission at these 
frequencies, an observed lower limit on R corre- 
sponds to a lower limit on /3 and an upper limit 
on i. The GBI radio light curves during the pe- 
riod when these observations showed a jet suggest 
a ~ —0.6, typical for optically-thin synchrotron 
emission; since the jet appears fairly continuous, 
wc take k = 2. i? > 330, as discussed above. We 
have then that 

/3cosi > 0.806 

implying /? > 0.81 and i < 36°. These limits arc 
consistent with, and more stringent than, those 
based on the proper motion. 

Clearly it would be useful to place limits on 
the maximum as well as the minimum proper mo- 
tion. Each epoch's observations took place over 
13 hours; one could therefore hope to check for 
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changes in source structure within a given observ- 
ing run. Unfortunately the hmited uv-coverage 
within a small segment makes it virtually impos- 
sible to image such a complex structure as the 
curved jet; further, the corresponding point spread 
function changes drastically as the earth rotates, 
making even an unchanging source appear quite 
differently when imaged with short segments of 
data separated by even a half hour. 

3.2.2. .Jet Curvature and Precessing Jet Models 

The simplest explanation for the curvature of 
the jet, and the change in that curvature between 
the first and the second image, is that the jet is 
precessing (cf. Hjellming & Johnston 1988). The 
fact that the jet curves out to the west, then back 
to the east, with no sign of another westward bend, 
implies that the ejection occured over at most one 
full precession period. Given the length of the 
jet this implies a joint limit on the speed and the 
period: 

P (cpsini) > I jet 

with Ijet being the length of the jet, 6.7 light-days 

on 8 February, assuming a distance of 10 kpc. The 
precession period P must therefore be 

P > Ijet/ {(Jsinic) « 6.7 days//? sin i 

Since /? sini < 1, P > 6.7 days. 

To proceed further requires more detailed mod- 
eling. Given the complex radio light curve and the 
various imaging problems mentioned above, and to 
avoid the additional assumptions involved in phys- 
ically modeling the ejecta, we attempt to fit only 
the shape of the jet and its evolution between the 
two epochs, not the jet's brightness distribution. 
That is, we only require that the models trace 
the jet geometry in both the first and the second 
epoch. We further assume, as in the last para- 
graph, that the observed jet was produced within 
a single precession period. 

The precessing jet model as described by 
Hjellming & Johnston (1988) has eight impor- 
tant parameters: the distance, d; the precession 
period, P; the intrinsic velocity, /3c; the inclina- 
tion angle of the jet axis to the observer, ija, the 
cone opening angle, ip] the position angle of the 
jet on the sky, P. A.; the phase of the precession 
at some fiducial time, (j^o; and the sense of ro- 
tation of the jet (clockwise or counter-clockwise. 



as seen from the jet origin). The distance of the 
source is kept fixed at 10 kpc; the other parame- 
ters are to be found by fitting the observed VLBI 
images. In order to find the full range of allowed 
parameters, we searched for the best alignment 
between the images and the model, based on 
minimization (Press et al. 1986), using several 
million randomly-chosen sets of input parame- 
ters. The initial guesses uniformly covered the full 
range of physically meaningful model parameters 
(see Table 2, with the maximum period (600 days) 
chosen as an approximation of 'very long'. The 
agreement between the model and the data was 
measured by the weighted sum of the squares of 
the distances between the closest points in the 
model, and ^ 10 positions measured along the 
'spine' of the jet in each of the first and second 
epoch images (see Figure 5a). The weights were 
taken from rough error bars based on the believ- 
ability of a feature and the approximate accuracy 
of its position, taking into account uncertainties 
due to noise in the images, deconvolution artifacts, 
the local width of the jet, etc. While the absolute 
value of the resulting is not very meaningful, 
models with lower do match the images better, 
which is all that is necessary for the minimization. 

The initial parameters, and the positions and 
error bars measured off the two images, were 
passed to a x^-minimization routine based on the 
downhill simplex method (Press et al. 1986), pro- 
ducing about a million converged solutions. Based 
on those results several hundred thousand more 
models were run, using a more restricted set of 
inputs chosen to maximize the number of 'good' 
solutions (defined below) , in order to flesh out the 
range of acceptable model parameters. We have 
three constraints on these solutions. First, they 
must match the observed jet morphology. Empir- 
ically x^ ^ 8 corresponds to a good chi-by-eye fit. 
Second, if the extended jet did indeed originate at 
MJD 50482. libO.l, the jet as seen in the first im- 
age can be at most 4.0 days old (as seen by us), 
while the jet seen in the second image can be at 
most 6.0 days old. Finally, if the jet is intrinsically 
two-sided, the jet/counter-jet ratio on 8 February 
is > 330 (see §3.2.1). Without a model for the 
jet brightness distribution and its evolution, this 
last constraint is impossible to apply fully to these 
models. We take a conservative approach, requir- 
ing that at least one observed point along the jet 
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Fig. 5. — Images of 6 February and 8 February observations with, contour levels of 0.025 x (1, 1.4) mJy/beam 
for the 6 February image and 0.0012 x (1, 1.4) mJy/beam for the 8 February imago, and a. the positions 
of the points used for in the model fits superimposed. The size of the crosses are 3cr and show the relative 
weight of each point (i.e., the smaller crosses have more weight in the model fit), b. a jet model with 
parameters of P = 29.8 days, (3 = 0.990, ija = 82? 1, tp = 55?7, counter-clockwise rotation; = 6.95 
superimposed. This model corresponds to the high-inclination branch of fits satisfying the age constraint 
(MJDgj > 50482.0). The corresponding peak Doppler ratio R is only 2.3 however, so these fits would 
require an intrinsic jet asymmetry or obscuration to explain the lack of a counter-jet. c. a jet model with 
parameters of P = 116 days, f3 = 0.902, ija = 6?3, ip = 1?1, clockwise rotation; = 7.99 superimposed. 
A barely-acceptable fit, this is an example of the clockwise-rotating family of solutions which satisfies the 
boosting but not the age constraint (it requires a minimum jet age of 7.1 days on 8 February), d. a 
jet model with parameters of P = 496 days, /? = 0.967, ija = l'?5, V = 1-2, counter-clockwise rotation; 
X^ = 6.44 superimposed. The best-fitting model of the family which satisfies all three (x^, boosting, and 
age) constraints. 



have a boosting factor, relative to the correspond- the images with examples of the major families of 

ing component in the purported counter-jet, of at solutions superimposed. 

least 3301/2-6 = 9,3 (gee §3.2.1). Figure 5b-c shows The results are given in Table 2. Good fits, as 
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measured by : may be obtained for a wide range 
of parameters (see Figures 6 & 7). In particular, 
jets pointing away from the observer can match 
the observed morphologies in both images. There 
is a lower limit on the period of about 15 days, 
analogous to the analytic limit discussed above. 
More surprisingly, there is also a lower limit on 
the age of the ejecta: they must have been ex- 
pelled at least 2.5 days before the first image, and 
at least 3.5 days before the second. Note that 
this is derived simply by fitting the morphology^ 
without any additional constraints from the radio 
light curves, and results from the requirement that 
both images be fit simultaneously within one pe- 
riod, which limits the proper motion and therefore 
the age. Interestingly, this lower bound on the age 
corresponds to the peak in the radio light curves 
around MJD 50483.5. 

Requiring in addition that the jet in the first 
image be no more than 4.0 (and the second no 
more than 6.0) days old, restricts the possible so- 
lutions enormously. Jets pointing away from the 
observer are entirely eliminated, because they can- 
not produce a long enough jet in the required time. 
Solutions with clockwise rotation (as seen from the 
jet origin) are also excluded, as they give poor fits 
to the observed morphology (minimum = 8.5). 
The solutions in which the jet rotates counter- 
clockwise are themselves split into two basic fam- 
ilies: low-inclination jets {ija ^ 20°) with a linear 
relationship between inclination and cone opening 
angle, and relatively long periods (P > 50 days); 
and high-inclination jets {ija > 50°) with large 
cone opening angles (35 — 90°) and relatively short 
periods. The latter branch is eliminated if suffi- 
cient boosting is required to conceal the counter- 
jet in an intrinsically symmetric system. The so- 
lutions obtained using all three constraints then 
form a one-parameter family, with 

rotation counter — clockwise 

P > 60 days 

/3 « (-11. 89/P) + 0.989 >0.81 

ija « (785?5/P) - 0.25° < 14° 

Tp w 0.79 - 0.08° < 12° 

P.A. « 145 ± 5° 

(ho « 70 ± 10° 




Fig. 6. — Inclination of the jet axis, ija, vs. cone 
opening angle, -0, for 'good' model fits to the 6 and 
8 February 1997 images. For simplicity, only those 
fits with counter-clockwise jet rotation and i < 90° 
are shown. Black dots represent fits with < 8. 
Requiring in addition that jet ejection begin no 
earlier than MJD 50482.0 (see text) restricts pos- 
sible fits to those indicated by the green (and blue) 
circles. Similarly, requiring that the jet's appar- 
ent asymmetry be due to Doppler boosting re- 
sults in the possible fits shown by the red (and 
blue) circles. Imposing both age & asymmetry 
conditions at once gives the single-parameter fam- 
ily of solutions shown as the blue circles. A least- 
squares fits to these latter points gives the blue 
line: ^/^ w 0.79 i^a - 0.08°. 



With these additional constraints, the minimum 
age for the first image is 2.8 days, corresponding 
to ejection on MJD 50482.7-50483.2. The limit on 
f3 arises directly from the jet/counter-jet ratio; the 
maximum ages for the two images constrain Papp, 
which is a combination of (3 and i; those in turn 
constrain the period, assuming the jet seen in the 
image was formed within a single precession pe- 
riod. The requirement to match the observed cur- 
vatures reduces this two-dimensional space of so- 
lutions to a single dimension, and further requires 
the above restrictions on the cone opening angle 
ijj and the other geometric properties {P-A., (f>o). 
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Table 2 

Results of Modeling Jet Morphology 





Rotation 


Period 
(days) 








P.A. 


00 


Initial guess'' 


CW, CCW 


- eoo'^ 


0-1 


- 180° 


0-90° 


- 360° 


- 360° 


<8 




>io 












< 8 & MJDej > 50482.0'= 


CCW 
CCW 


> 50 
20-60 


> 0.75 

> 0.75 


< 18° 
50 - 85° 


< 15° 
35 - 90° 


120 - 160° 
210 - 275° 


50 - 100° 
265 - 325° 


< 8 & > 330 


CCW 
CW 


> 60 

> 70 


> 0.81 

> 0.81 


< 13° 

< 15° 


< 18° 

< 10° 


60 - 300° 
180 - 220° 


> 50° 
135 - 270° 


< 8, MJDej > 50482.0? 
<t /? > :',:!() 


CCW 


> GO 


> ().<Si 






125 - i.-).-)'-' 


")0 - 



^Inclinations larger than 90° correspond to one-sided jets pointing away from the observer. 

''Initial guesses were taken as uniform random deviates over the listed range, with each parameter chosen indepen- 
dently. 

°Both clockwise (as seen for the jet origin) and counter-clockwise jet rotation were allowed. 
"^The maximum period (600 days) was chosen as an approximation of "very long." 
''MJDej is the date at which jet ejection began. 
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0,03 0,04 

1/P [days-i] 







0,03 0,04 
1/P [days-i] 



Fig. 7. — Inclination of the jet axis, ija (left), and intrinsic jet speed, P (right), as a function of the 
inverse period, 1/P. Colors are as in Figure 6. The blue lines indicate least-squares fits to the blue circles: 
ija « (785?5/P) - 0.25°; /3 « (-11.89/P) + 0.989. The gap on the left-hand side of each plot corresponds to 
the maximum period allowed as an initial guess, 600 days. 



In sum, the simple precessing-jet model of 
Hjellming & Johnston (1988) can fit the observed 
jet quite well, for either a one- or a two-sided jet. 
To conceal the counter-jet requires substantial 
Doppler boosting — in these models the intrinsic 
(as compared to the differential jet/counter-jet) 
Doppler boosting of the flux is at least a factor 
13, and could be much higher. Also, for these low- 
inclination models, temporal variations in the jet's 
frame are significantly compressed in ours because 
of time dilation (e.g., in one model 6 months of jet 
evolution looks like only 6 days to the observer of 
the jet). This might help explain both the unusual 
strength and the rapid variability of Cygnus X-3 
as a radio source, although it should be noted that 
this time dilation would not affect variability that 
originates in the accretion disk or any other part of 
the system that is not moving towards us at rela- 
tivistic speeds. The small range of position angles 
observed over the years is also consistent with the 
small range of ejection angles in these models - jet 
models more nearly in the plane of the sky require 
much larger cone angles to fit our images, giving 
correspondingly wider position angle swings over 
the full precession period. 



One argument against the low-inclination mod- 
els is that the periodic dips in the X-ray and in- 
frared light-curves are usually modeled as opacity 
effects in a reasonably edge-on disk (e.g., Mason, 
Cordova, & White 1986). If the radio jet were per- 
pendicular to the disk this would suggest the jet 
must be in the plane of the sky. However, Ghosh 
et al. (1981) showed that the X-ray light curves 
could be matched in detail even for disks with 
inclinations up to 70°, and observations at other 
wavelengths are also consistent with a wide range 
of inclinations (Schmutz, Geballe, & Schild 1996; 
Hanson, Still, & Fender 2000). Also, we clearly 
see precession in the jet, which implies that the jet 
is not perpendicular to the binary system's plane 
of orbital motion. So our solutions do not put a 
tight constraint on the inclination of the binary 
system. 

It is also possible that the system is edge-on and 
the jet in Cygnus X-3 is in the plane of the sky. 
The one-sided jet could then be explained either as 
intrinsic or, as suggested by Fender et al. (1999b), 
by obscuration. Obscuration seems very unlikely, 
since it would require not only an odd geometry, 
blocking out one side of the jet but not the core or 
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the other side of the jet, but also material opaque 
at 15 GHz out to > 1200 AU. 

4. Conclusions and Implications 

The main result of these observations is to sliow 
that the radio emission from Cygnus X-3 during 
flares is dominated by a one-sided relativistic jet 
with an intrinsic speed of at least 0.81c. Assuming 
the jet to be intrinsically symmetric, precessing jet 
models give a maximum inclination to the line-of- 
sight of ~ 14°. Cygnus X-3 is the most luminous 
X-ray binary at radio wavelengths, and the ob- 
served asymmetry and high proper motion make 
it tempting to speculate, as mentioned in the pre- 
vious section, that this high luminosity is due to 
significant Doppler boosting of a jet pointed al- 
most directly towards us. This might also help 
explain the rapid and extreme radio variability of 
this and other known relativistic jet sources. How- 
ever, since no counterjet is seen, the VLBA ob- 
servations do not require such boosting, and, un- 
less the X-rays are also associated with the jet 
(contrary to most current models), boosting of 
the jet would not explain why Cygnus X-3 is also 
one of the brightest X-ray sources in the Galaxy 
(10^'^ - 10^^ erg/s, Bonnet-Bidaud & van der Klis 
1981). 

There are three other relativistic jet sources in 
the Galaxy which have been studied in some de- 
tail: SS433, GRS 1915+105 and GRO J1655-40. 
Cygnus X-3 is the only one-sided jet among the 
four, presumably because at the time of the flare 
it was aligned more nearly along the line-of-sight, 
although the jet GRO J1655-40 is sometimes in- 
trinsically asymmetric (Hjellming & Rupen 1995). 
Like GRO J 1655-40, Cygnus X-3 is a strong black 
hole candidate based on its mass function, a con- 
clusion strengthened by the limit derived here 
on the inclination of the jet, which is probably 
roughly aligned with the angular momentum axis 
of the orbit. Cygnus X-3 is however the flrst con- 
sistently strong X-ray source to exhibit such highly 
relativistic jets; GRS 1915+105 and GRO J1655- 
40 are both X-ray transients, often undetectable 
but occasionally among the brightest sources in 
the sky, while SS433 is at best an undistinguished 
X-ray source. Similarly, only Cygnus X-3 and 
SS433 are detectable in the radio even wIkui they 
are not flaring, though GRS 1915+105 can remain 



in a plateau state for months or longer (Foster 
et al. 1996). The radio and X-ray states are 
closely coupled for all three of the highly relativis- 
tic jets, while for SS433 the lack of an obvious con- 
nection may simply be due to the relatively poor 
X-ray coverage. Intriguingly, both GRO J1655- 
40 and GRS 1915+105 have a very imusual, hard 
power-law X-ray tail to energies of several hun- 
dred keV, while Cygnus X-3 has become famous 
as an occasional source of these (Matz et al. 1996) 
and even higher energy photons (certainly up to 
100 MeV; see Lamb et al. 1977 and Fichtel, 
Thompson, & Lamb 1987). 

The clear implication is that changes in the ac- 
cretion disk produce changes in the radio jet, and 
that an excess of high-energy photons may indi- 
cate a source capable of producing highly relativis- 
tic radio jets. It is not particularly surprising that 
the high-energy photons and the high-energy elec- 
trons should be fairly closely coupled. What is 
new here is that an X-ray binary, with consistent 
and reasonably strong X-ray and radio emission — 
and hence, a relatively stable accretion rate, and 
accretion disk — can give rise to the highly rel- 
ativistic jets previously associated only with very 
unusual X-ray transients. 

An obvious question is whether the quiescent 
radio emission is also in the form of a relativistic 
jet. In SS433, it clearly is, and neither the intrin- 
sic speed nor the orientation of the jet depend on 
the strength of the X-ray or radio emission. Nei- 
ther of the highly relativistic X-ray transients has 
been detected in radio quiescence, although obser- 
vations during the smaller flares of GRS 1915+105 
are broadly consistent with the long-lived flare im- 
aged by Mirabel & Rodriguez (1994). Unfortu- 
nately the current data on Cygnus X-3 are still too 
fragmentary to convince one way or the other. Our 
11 February observations suggest at most a much 
slower expansion rate than during giant flares, 
while NGS find something much faster. Despite 
this confusion Cygnus X-3 offers the unique oppor- 
tunity to check the behavior of a highly relativistic 
jet source over its whole range of X-ray and radio 
states; we and doubtless others will be observing 
this source for some time to come. 

Finally, with a convincing jet found in an X-ray 
binary like Cygnus X-3, it is beginning to seem 
that every Galactic X-ray source with radio emis- 
sion turns out, when imaged, to be a relativistic 
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jet. While the four discussed so far arc the only 
truly compelling examples, there are hints that 
several others are jets as well (possible polariza- 
tion in 4U 1630-47, Buxton ct al. 1998; elongated 
structure in GX 339-4, Fender et al. (1997); sug- 
gestions of elongation and high-speed expansion in 
LSI-H61°303, Pcracaula, Gabuzda, & Taylor 1998; 
VLA images of one-sided jet in V4641 Sgr, Hjellm- 
ing et al. 2000). On the other hand, all the rela- 
tivistic sources currently known are very peculiar 
in some ways, although few of those peculiarities 
are the same for all four sources. The next chal- 
lenge is to image some of the fainter, more common 
radio X-ray binaries, to see whether more 'normal' 
systems also give rise to relativistic jets. 

This project could not have been done without 
the Green Bank Interferometer (GBI), which pro- 
vided the flux measurements which triggered our 
target-of-opportunity observations. Those obser- 
vations themselves could not have succeeded with- 
out the eager help of many people associated with 
the VLBA, particularly the site techs who ran up 
various mountains to change tapes at very odd 
hours; we are most grateful both to them and to 
the scheduling committee for their timely efforts. 
M.R. enjoyed a month's hospitality at the RCfTA 
at Sydney University. Finally, Ketan Desai kindly 
provided, discussed, and improved his uv-plane fit- 
ting program in response to our requests. The GBI 
is a facility of the National Science Foundation op- 
erated by the National Radio Astronomy Obser- 
vatory (NRAO) in support of NASA High Energy 
Astrophysics programs. NRAO is a facility of the 
National Science Foundation operated under coop- 
erative agreement by Associated Universities, Inc. 
A.M. acknowledges support from the European 
Commission's TMR/LSF Programme (Contract 
No. ERB-FMGE-CT95-0012). Basic research in 
radio astronomy at the Naval Research Labora- 
tory is funded by the Office of Naval Research. 
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